Disorder-induced zero-energy spectral singularity for random matrices with 

correlations 
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A zero-energy mid-band singularity has been found in the energy spectrum of random matrices 
with correlations between diagonal and off-diagonal elements typical of vibrational problems. Two 
representative classes of matrices, characterizing the instantaneous configurations in liquids and 
mechanically unstable lattices (which mimic the former) have been analysed. At least for disordered 
lattice models, the singularity is universal and its origin can be explained within the mean-field 
treatment. 



Many physical phenomena can be described by the 
Anderson Hamiltonian with ofF-diagonal disorder (see 
e.g. and references therein). Stochastic transport 0, 
atomic vibrations in disordered structures (see e.g. B] 
and references therein) and instantaneous normal mode 
(INM) analysis in liquids and glasses Q| (called below "vi- 
brational" problems) arc among these. However, there is 
one key point which distinguishes vibrational problems 
from the standard electron problems: there exist strong 
sum-rule correlations between the ofF-diagonal and diag- 
onal elements (the sum of all the elements for a particular 
row is zero) of the relevant random matrices for the vi- 
brational problem. 

In different dimensions (D = 1 — 3), the energy spec- 
trum of the electron Anderson Hamiltonian with pure 
ofF-diagonal disorder defined on the simple (hyper) cubic 
lattice (but not, as we have checked, on the f.c.c. lattice) 
exhibits a peculiar feature: a disorder-induced, mid-band 
(zero-energy) singularity occurs I, §. Normally, disor- 
der smears out sharp features in a spectrum (e.g. van 
Hove singularities), but in this case it creates a singu- 
larity, the origin of which in different dimensions is still 
controversial. However, can disorder induce a similar (or 
different) singularity for vibrational problems and, if so, 
can its origin be revealed? In this Letter, we demon- 
strate the existence of a mid-band zero-energy singular- 
ity in the spectrum of dynamical matrices involved in an 
INM analysis of a model of a monatomic liquid, and in- 
vestigate analytically and numerically the nature of this 
zero-energy singularity in the spectrum of random matri- 
ces with sum-rule correlations, appearing in disordered 
lattices which mimic topologically disordered liquids. 

We have observed the occurrence of a singularity by 
numerical experiment, in which structural models of a 
liquid with predominantly icosahedral order (similar to 
those discussed in Ref. Q) and of glassy and liquid silica 
(similar to those discussed in Ref. ||l^) have been con- 
structed at different temperatures, T, by slow molecular 
dynamics quenches (by constant volume and tempera- 
ture steps). The INM spectra of the dynamical matrix 
are shown in Fig. A zero-energy singularity is clearly 
present in the energy spectrum, g{£), for the icosahedral 
liquid and possibly in silica (at least, the tendency of the 
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FIG. 1: The energy spectrum, g{e), for the instantaneous 
dynamical matrices for an INM analysis of (a) a liquid with 
predominantly icosahedral order (1620-particle model aver- 
aged over 90 configurations) at two temperatures, as marked 
(the units are the same as in Ref. [^]) and of (b) glassy (the 
solid curve) and liquid (the dashed curve) silica (1650-particle 
model averaged over 20 configurations). The energy units for 
silica are THz^. 



change of the DOS around zero is consistent with the 
model presented below) . Similar zero-energy peaks have 
been found hut not explained in liquid theories 10,^ and 
in simulations of a binary Lennard- Jones liquid . Nor- 
mally, the density of states (DOS) in the frequency do- 
main, — 2yJ\£\g{£), is of interest in the INM analysis, 
but the singularity is masked there by the linear factor 
2w, and that is, probably, why it has not been carefully 
studied before. 

Therefore, there exists a peculiarity in the spectrum 
of random matrices with sum-rule correlations (instanta- 
neous dynamical matrices) and its origin can be revealed 
analytically as follows. Consider a Hamiltonian describ- 
ing both the electron and vibrational problems: 



H 



{e^^lY.Uj)\i){i\^Y. ^-o\i){3\ , (1) 



where Si and tij stand for the random on-site energies 
and the random transfer integrals between sites i and j, 
respectively. The parameter 7 controls the correlations 
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FIG. 2: (a) Evolution of the zero-energy singularity in the 
DOS for a disordered f.c.c. lattice with various values of dis- 
order, A (A* ~ 1.3 and A«« ~ 2.2). The thick curves were 
obtained by the numerical integration of the CPA DOS given 
by Eq. (^), while the thin ones are the results of precise nu- 
merical KPM solution for vector vibrations in a 130 x 130 x 130 
site f.c.c. lattice. The difference between the CPA and KPM 
curves around both band edges is due to the known failure of 
CPA to reproduce the localized band tails. The inset shows 
the phase, <?!>(0), of the effective field at zero energy versus 
disorder, A. The critical values of disorder are indicated by 
arrows, (b) An enlargement of the zero-energv region. The 
CPA results obtained by use of Eqs. and are not dis- 
tinguishable on this energy scale. 

between diagonal and off-diagonal matrix elements. The 
standard electron Anderson Hamiltonian with pure off- 
diagonal disorder corresponds to 7 = and Ei — 0, with 
tij being random variables taken, for example, from a 
uniform (box) distribution of width 2A centred around 
to = —1, tij S [to — A, to + A]. The vibrational (scalar) 
problem corresponds to 7 = 1 with the other parame- 
ters being the same. For atomic vibrations, the Hamil- 
tonian corresponds to the dynamical operator, transfer 
integrals to the force constants and energy to the squared 
frequency, e — lu'^ 

We have calculated by the kernel polynomial method 
(KPM) the spectrum of the "vibrational" Hamilto- 
nian (IJ) (see the thin lines in Fig. ^) for 7=1 and 
Ei = 0, defined on the f.c.c. lattice (actually, of its more 
general version for vector vibrations; see |^ for more de- 
tail). For sufficiently large degrees of disorder, A > Aq 
(Aq ~ 1.4 — 1.5), the zero-energy singularity is evident 
(Fig. H). Moreover, the shape of the singularity is simi- 
lar to that found for the topologically disordered models 
(Fig. ^, indicating its possible universality (see below). 
We have also calculated the DOS for the same problem 
within a mean-field approach (the single-bond coherent 
potential approximation (CPA) |^), and found remark- 
ably good agreement with the precise numerical (KPM) 
results (cf. thick and thin lines in Fig. ^ in the singular- 
ity region for A > Aq, ((Ao — A*)/A* ^ 1, see below for 



a discussion of A*), i.e. when the localization threshold 
is far enough below zero energy. This agreement is sur- 
prising (at first sight), because it is commonly believed 
(see e.g. [Q) that mean-field theories fail to reproduce 
sharp features in a spectrum. For example, the most suc- 
cessful homomorphic cluster CPA |l^ well reproduces 
the whole spectrum for the electron problem with pure 
off-diagonal disorder (7 = 0), except for the zero-energy 
singularity. Below, we explain why the CPA reproduces 
the zero-energy singularity for vibrational problems but 
not for the electron one, and we use this insight to reveal 
the physical origin of the singularity. 

The zero-energy singularity occurs for 7=1 and disap- 
pears if 7 7^ 1, i.e. when the exact sum-rule correlations 
of elements in the dynamical matrix break down. The 
functional form of the zero-energy singularity observed 
is universal (i.e. independent of the type of distribution 
of the force constants, the reference lattice symmetry, 
scalar or vector type of vibrations, etc.) for the class of 
Hamiltonians given by Eq. (|l|) with 7 = 1 and Si = 0, 
and depends only on the dimensionality of the problem. 
For sufficiently large disorder, the (parity-independent) 
singularity occurs in the mid-band region of the spec- 
trum, far from the localized states in the 3Z?-case (as 
checked by multifractal analysis) and in the range of pre- 
localized states for lower dimensions. The analytic mean- 
field solution shows that its appearance is dictated by the 
universal non-analytic behaviour of the spectral-density 
operator in the plane- wave basis and, essentially, is a con- 
sequence of the fact that zero-energy (long-wavelength) 
plane waves contribute anomalously (but not solely) to 
the disordered eigenstates with energies around zero. 

The DOS, g{e), for a disordered lattice is the trace 
of the spectral-density operator, A(e) = {S{e — H)), 
taken in the convenient orthonormal basis of crys- 
talline eigenstates |k, /3) (k is the wavevector of a plane 
wave from the branch /3) of the same Hamiltonian, 
but without disorder (A — 0), and averaged ((. . .)) 
over the distribution of random variables tij, g{e) — 

/((k,/3|A(£)|k,/?))5-y^t(£k/3)dek^, with ^"'''^(ek/?) be- 
ing the crystalline VDOS. The diagonal matrix el- 
ement of the spectral-density operator, yl(e,£k/3) = 
(k,/3|A(e)|k,/3), can be found within the CPA |, |l|, 
[Tsl , [T^ via the complex effective interaction field z(e) — 
z'{e) + \z"{e), so that 

n(r\ - ^ [ (gk/3/£)r(£) „cryst/^ N 

(2) 

with 

e(e) = ez'(e)/|z(£)p and r(e) = ~ez"ie)/\z{e)\^ . (3) 

The effective interaction field, ^(e), is found from the so- 
lution of the self-consistent CPA equation (see |^). The 
real part, z'{e), fluctuates around its crystalline value 
Zq = 1, and the imaginary part is smooth and negative 
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FIG. 3: Dependence of the spectral density, (j4(e, Sk/j)), on 
the crystaUine energy, £k/3, at values of disordered energy, e, 
as marked (for A = A**), for vector vibrations in the force- 
constant disordered f.c.c. lattice. The inset shows the en- 
ergy dependence of the spectral-density peak width parame- 
ter, r(£) (Eq. (0)), for different values of disorder, as marked. 



in the region of the non-vanishing disordered spectrum, 
and zero otherwise. The spectral density, (A(£,£k/3)}, as 
a function of ek/3, has the shape of a peak, character- 
ized by a width parameter, r(e), located at an energy, 
Epcak = |e|k(e)r' = [e{e)+T{e)]'^^ (see Fig.|). 

As follows from Eq. (||), the disordered DOS is gener- 
ically related to the crystalline DOS for the reference 
system convolved with the peak-shaped spectral density, 
{A{e, etifj)) . For small disorder (A 0), the spectral 
density is very narrow and close in functional form to 
S{£ — £kp), SO that the disordered DOS strongly resem- 
bles its crystalline counterpart. With increasing disorder 
but for A < A* (where the spectrum of the Hamilto- 
nian is still non-negative and the lattice is mechanically 
stable ^), the spectral-density peak becomes broader 
and washes out all the van Hove singularities in the crys- 
talline spectrum, except the boundary singularity around 
the lowest band edge, Smin = 0. 

When the disorder exceeds the critical CPA value, 
A > A*, the number of negative force constants becomes 
so large that the system is no longer mechanically stable 
even in the mean- field limit, and negative eigenvalues 
appear in the spectrum (as in instantaneous configura- 
tions of liquids Q), i.e. the lower boundary of the CPA 
spectrum moves below zero energy, £niin < 0. Exactly 
in this regime, the mid-band zero-energy singularity in 
the disordered CPA spectrum evolves. Indeed, if we look 
at the evolution of the spectral-density width parameter, 
r(£), with increasing disorder (see the inset in Fig. ||), 
we can clearly see that the peak width is still zero at 
zero energy, r(0) = 0, being a consequence of Eq. (||), 



even though the value of the effective field becomes finite 
at this point {z"{Q) ^ 0). This means that the spectral 
density, {A{e — 0, £k;3)), as a function of the crystalline 
energy, £k/3, has a singularity at £k/3 = (see Fig. ||). 
The zero-energy point, in this regime (A > A*), belongs 
to the mid-band region. The finite value of z"(0) im- 
mediately gives rise to a different singular shape of the 
spectral density, (A(0,£k/3)) ~ -7r-i(z"(0)/|z(0)|2)£-^ 
(see Fig. ||), as compared to the 5- functional shape of 
the spectral density for the A < A, regime, where both 
r(0) and z"(0) are zero. The (S-functional shape of the 
spectral density reproduces the crystalline van Hove sin- 
gularity, g(e) oc £(^/^)~^, at the lower band edge, but 
the £^^-shape of the spectral density at £ = 0, convolved 
with the crystalline DOS (see Eq. (^), results in a sin- 
gularity of a different type in the disordered DOS which 
is not related to the crystalline van Hove singularity at 
£ = (see below). 

The shape of the mid-band zero-energy singularity can 
be obtained analytically within the CPA by splitting the 
integration region in Eq. (^ into a low-energy part, £k/3 < 
£o, where g'^'^^^^{e) ~ Xd£^^^^^~^ (the Debye law), and 
the rest of the band, £k/3 > £o (irrelevant for the shape of 
the singularity). The final result for the disordered DOS 
(3D-case) in the limit |£| ^ is the following: 

Csinc/)(£) , xM^I^ . 

(4) 

where C = J [g"y^^ (£) /e] de is a model-dependent con- 
stant, (j){e) = arg [z{e)] is the phase of the effective field 
and 9{e) is the Heavyside function, 9{e < 0) = and 
6'(£ > 0) = 1. Bearing in mind that the effective field is a 
smooth function of energy around zero, we can conclude 
from Eq. (^) that, in the 3D-case, the disordered DOS 
is a continuous function around zero but its derivative 
shows £~^/^ singular behaviour (see Figs. 0, ||). 

The expression (||) for g{e) has been obtained in the 
limit |£| for any degree of disorder. The evolution of 
the singularity with disorder (see Fig. ||) is determined 
by the dependence on A of the effective field at zero en- 
ergy. The phase, 0(0; A), varies in a critical manner with 
disorder (see the inset in Fig. §(a)), being zero below the 
critical disorder, A < A^ , when the zero-energy singular- 
ity is just a lower band-edge van Hove singularity. Above 
the critical disorder, A > A*, the derivative of the DOS 
becomes singular from both sides about £ = 0. If the 
disorder is not too large. A* < A < A*,, the sign of the 
derivative is positive on both sides of the singularity and 
the DOS is monotonic in the singularity range (see the 
upper curve in Fig. §(b)). For higher disorder, A > A*,, 
the derivative changes sign at the singularity and the dis- 
ordered DOS exhibits a sharp maximum at zero energy 
(see the lower curve in Fig. ^(b)). The characteristic 
value of disorder, A = A*,, at which such a transforma- 
tion of the shape of the singularity occurs (the solid line 



4 



in Fig. g), can be found from the solution of the equation 
(/)(0;A**) = — 7r/3 (see the inset in Fig. ||(a)), resuhing 
from the condition, g'{e 0+) = 0. 

The physical significance of this transition is related to 
the fact that, at A ~ A*^,, the peak width of the spec- 
tral density around the singularity becomes comparable 
with the peak position. This means that A*, corresponds 
to the loffe-Regel crossover for the propagation of plane 
waves characterized by the energies e^p In 
other words, for A < A*^,, there is a finite low-energy in- 
terval of the weak-scattering regime for plane- wave prop- 
agation, but for A > A»*, all plane waves propagate in 
the regime of strong scattering. 

Similar CPA analyses can be performed for lower di- 
mensions, and they result in zero-energy singularities 
as well. In 2D (for A > A,), a logarithmic singu- 



larity, g{\£\ 0) 



^(2"(e)/|z(£)|2)x2ln|£|, evolves 



on the background of the van Hove band-edge singu- 
larity. In ID, the divergence of the disordered DOS 
from both sides of zero energy is even more pronounced: 
g{e) c ^xi\e\-'^Me)\-'^'Mm/2-0{e)7r/2). 

The link between lattice models and topologically dis- 
ordered liquids and glasses is, of course, not straight- 
forward. The above mean-field analysis for disordered 
lattices uses the existence of an orthonormal plane-wave 
basis, which is absent in topologically disordered systems. 
Nevertheless, an orthogonalized basis resembling a plane- 
wave basis, at least in the region of zero energy, can read- 
ily be constructed , so that it might be expected that 
the same universal singular behaviour of the DOS around 
e = should also occur for topologically disordered sys- 
tems. Indeed, the singularity is evident for a model liquid 
with predominantly icosahedral order (Fig. |l|(a)), and it 
is remarkable that the evolution of the shape of this sin- 
gularity with decreasing temperature mirrors that found 
in the lattice models with decreasing force-constant dis- 
order (Fig. |^(b)). Both regimes of the singular behaviour 
are evident from the numerical experiment (Fig. 0(a)): (i) 
the high-temperature liquid state (circles) can be char- 
acterized by relatively large disorder (A > A**); (ii) the 
state just above the glass-transition temperature (solid 
curve) can be associated with the intermediate disorder 
regime (A* ^ A < A*»). The situation for vitreous 
silica (Fig. 0(b)) is different. At both temperatures, be- 
low and above glass transition, the models stay in the 
regime A < A*,, which is not surprising because silica 
is a very strong glass-forming liquid (especially in com- 
parison with the fragile system shown in Fig. 0(a)), and 
is characterized by well-defined and stable local tetrahe- 
dral order even at very high temperatures (this means 
that the disorder, i.e. A, is relatively small). A more de- 
tailed analysis (including the glass-transition region) will 
be presented elsewhere. 

In conclusion, we have demonstrated the presence of a 
mid-band zero-energy singularity in the spectrum of the 
instantaneous dynamical matrices of topologically disor- 



dered structural models of liquids and of the dynamical 
matrices of disordered lattices which well mimic the for- 
mer. The presence of exact sum-rule correlations between 
the diagonal and off-diagonal elements in the disordered 
dynamical matrix causes this singularity. The shape of 
the disorder-induced singularity, at least in lattice mod- 
els, is universal and depends only on the dimensionality of 
the model. Such a universality is related to the universal 
ej^^ -behaviour of the spectral density in the plane- wave 
basis, which is due to the multiplicative nature of the 
effective interaction mean field (the mean-field energy of 
the quasi-particles, ek/3, is the product of the dimension- 
less effective interaction, z{e), and the bare crystalline 
energy, i-e ek/3 = z{s)ekf3). This property distin- 
guishes the above class of Hamiltonians from those with 
pure on-site energy disorder, characterized by an additive 
effective field 10 , which do not exhibit a zero-energy sin- 
gularity. Thus, the failure of the homomorphic effective 
field (containing both additive and multiplicative contri- 
butions) in reproducing the zero-energy singularity for 
the Anderson Hamiltonian with pure off-diagonal disor- 
der without correlations in simple cubic lattices can be 
explained. A successful effective field for such a prob- 
lem should not contain the additive part, at least at zero 
energy, but its construction is still an open question. 
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